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Abstract 

t-h ■ Intermediate energy Coulomb excitation and dissociation is a useful tool for nuclear structure 

. and astrophysics studies. Low-lying strength in nuclei far from stability was discovered by this 

00 | method. The effective range theory for low-lying strength in one-neutron halo nuclei is summarized 

and extended to two-neutron halo nuclei. This is of special interest in view of recent rather accurate 
experimental results on the low- lying electric dipole strength in 11 Li. Another indirect approach 

\Q ', to nuclear astrophysics is the Trojan horse method. It is pointed out that it is a suitable tool to 

investigate subthreshold resonances. 



a 



"o 1 Introduction and Overview 

With the exotic beam facilities all over the world - and more are to come - direct reaction theories are 
experiencing a renaissance. These direct reaction theories are essential to interprete indirect methods 
^ ■ used in nuclear astrophysics. We discuss Coulomb dissociation ;lj and the Trojan horse method, for 
a recent review see Coulomb dissociation was reviewed in 2003 [3j, a comprehensive paper on the 
Trojan horse method appeared in A more recent (mini)review of indirect methods can be found in 
j^]. Here we report on the effective range theory of low-lying electromagnetic strength in one-nucleon 
halo nuclei jfllll], and its extension to two-neutron halo nuclei. This will be of special interest in view 
of recent rather accurate experimental results on the low-lying electric dipole strength in n Li jH]. 

The Trojan-Horse Method 4j is a particular case of transfer reactions to the continuum under 
quasi-free scattering conditions. Special attention is paid to the transition from reactions to bound and 
unbound states and the role of subthreshold resonances [SL Since the binding energies of nuclei close 
to the drip line tend to be small, this is expected to be an important general feature in exotic nuclei. 



2 Effective Range Theory of Halo Nuclei 

2.1 Single particle halo nuclei 

At low energies the effect of the nuclear potential is conveniently described by the effective-range 
expansion JU]. An effective-range approach for the electromagnetic strength distribution in neutron 
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halo nuclei was introduced in O [Z| and applied to the single neutron halo nucleus 11 Be. Recently, the 
same method was applied to the description of electromagnetic dipole strength in 23 [TTJ. A systematic 
study sheds light on the sensitivity of the electromagnetic strength distribution to the interaction in 
the continuum. We expose the dependence on the binding energy of the nucleon and on the angular 
momentum quantum numbers. Our approach extends the familiar textbook case of the deuteron, that 
can be considered as the prime example of a halo nucleus, to arbitrary nucleon+core systems, see also 
[T2| . We also investigate in detail the square- well potential model. It has great merits: it can be solved 
analytically, it shows the main characteristic features and it leads to rather simple and transparent 
formulae. As far as we know, some of these formulae have not been published before. These explicit 
results can be compared to our general theory for low energies (effective-range approach) and also to 
more realistic Woods-Saxon models. Due to shape independence, the results of these various approaches 
will not differ for low energies. It will be interesting to delineate the range of validity of the simple 
models. 

Our effective-range approach is closely related to effective field theories that are nowadays used 
for the description of the nucleon-nucleon system and halo nuclei |T3]. The characteristic low-energy 
parameters are linked to QCD in systematic expansions. Similar methods are also used in the study 
of exotic atoms {tt~A, 7r + 7r~, n~p, . . . ) in terms of effective-range parameters. The close relation of 
effective field theory to the effective-range approach for hadronic atoms was discussed in Ref. |14j . 

In the ANC method ^3] it is assumed that the radial continuum wave function gi is a pure Coulomb 
wave function (a regular spherical Bessel function for neutrons) and that the dominant contribution to 
the radial integral 

POO 

Ir = Jo ^( r ) r %( r ) rfr W 

comes from the region r > R where R is the range of the nuclear interaction. In this case the S-factor 
is entirely determined by the asymptotic normalization constant (ANC) of the bound state wave 
function 

fk(r) ~ C k e-« r (2) 

where U is the orbital angular momentum. The ANC is determined by a suitable transfer reaction. In 
some cases this is a good approximation. In general there is an influence from the final state interaction 
on the final continuum wave function. In the outside region (r > R) it is given by 

g tf = F if cos 5 if + G if sin S if (3) 

where the regular and irregular Coulomb wave functions are denoted by Fi f and Gi f respectively. For 
the low energies relevant now, the phase shift 8i f is given by the low energy parameters of the effective 
range expansion. Usually this is one number, the scattering length. 

In Fig.^we show the application of the method to the electromagnetic dipole strength in 11 Be. The 
reduced transition probability was deduced from high-energy 11 Be Coulomb dissociation at GSI [16]. 
Using a cutoff radius of R — 2.78 fm and an inverse bound-state decay length ofq = 0.1486 fm _1 as input 
parameters we extract an ANC of C = 0.724(8) fm" 1 / 2 from the fit to the experimental data. The ANC 
can be converted to a spectroscopic factor of C 2 S = 0.704(15) that is consistent with results from other 
methods. In the lowest order of the effective-range expansion the phase shift in the partial wave with 
orbital angular momentum I and total angular momentum j is written as tan 5/ = — (xc , l , y) 2l+1 . The 

halo expansion parameter is denoted by 7 = qR = 0.4132 < 1 and x = k/q = ^ E / S n where the neutron 
separation energy is given by S n . The effective range term \rik 2 term can be neglected, since it leads to 
a contribution with an extra 7 2 factor which is small in the halo nucleus limit 7 — > (at least in the case 
where the scattering length and the effective range parameter are of natural order). The dimensionless 
parameter corresponds to the scattering length a{ = (cjR) 2l+1 . We obtain c^ 2 = —0.41(86,-20) 
and Ci" 2 = 2.77(13, —14) , corresponding to a 1 ^ 2 = 456(61, — 73)fm 3 . This unnaturally large value was 
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Figure 1: The reduced electric dipole transition probability is shown as a function of the excitation 
energy E* = E + S n . The theoretically calculated dipole strength is compared to experimental data 
extracted from the Coulomb dissociation of 11 Be |16j . 

found directly from the analysis of Coulomb dissociation data JH] • It finds a natural explanation in the 
existence of the \ -bound state close to the neutron breakup threshold in n Be. 

The connection of the scattering length ai and the bound state parameter q for I > is given by 
ai = g 2(2^)H(2/-i)!! m a sc i uare we U model with radius R. This is a generalization of the well-known 
relation = 1/q for I = 0. The pi/2 channel in n Be is an example for the influence of a halo state on the 
continuum. The binding energy of this state is given by 184 keV, which corresponds to q = 0.094 fm -1 . 
With R = 2.78 fm one has 7 2 = 0.068. For I = 1 one has a\ = = 210 fm 3 . Applications to low-lying 
single neutron halo strength in carbon isotopes are given in [T7] . 

Let us now mention very briefly the proton case [Zj. Only a numerical approach is possible now. 
Again one finds that the ANC is not the only factor which determines the S-factor: there is also an 
influence of the scattering length. In Fig. 2 the influence of the potential depth Vq on the strength 
distribution is shown. We refer to [7] for further discussion. 

2.2 Two neutron halo nuclei, application to n Li 

While single particle halos show rather simple features which can largely be explained in the frame of 
the one-body Schrodinger equation, two neutron (or two proton, or neutron-proton) halo nuclei show 
many more interesting aspects. 

A Borromean nucleus is a three-body system which falls apart when one of the particles is removed, 
n Li is an example. It can be considered as a bound system of two neutrons and a core ( 9 Li). The 
system is barely bound (Eb ~ 300keV) and the core can be assumed to be inert. The subsystems (n-n) 
and ( 9 Li-n) are unbound. The binding of the Borromean nucleus is due to the attractive n-n interaction 
which leads to the binding of the three-body system, see e.g. [T8] . 

Essentially due to the small binding energy the low-lying El strength shows simple universal features 
as we will show now. The system is probed by long wavelength photons and can be characterized by 
a few low energy constants. As we have seen in the previous subsection, the dipole matrixelement is 
dominated by contributions from distances larger than the range of the core-neutron interaction. In 
this region the wave function depends on the binding energy Ej,. In a mean field approach the binding 
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Figure 2: Reduced transition probablity dB(El)/dE for the breakup of 8 B into a proton 
and 7 Be as a function of the relative energy E of the final state for various depths Vq of the 
potential in the continuum. The proton separation energy in 8 B is given by 0.1375 MeV. 
For further explanations see [7]. 
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energy is shared equally between the two neutrons. Each neutron is bound by e = ^p- = E b /2. The 
slope of the wave function is characterized by q, the overall magnitude by an asymptotic normalization 
constant Cf 0z, e (V) —> Ci£~ qr - 

For simplicity let us neglect spin-orbit forces and assume that the neutron spin state is the anti- 
symmetric S = state. The space part is symmetric , we include s 2 and p 2 configurations, with the 
corresponding spectroscopic factors C 2 Si. The radial part of the two-neutron halo wave function is 
written as ~ bi < t>i,e{ r i) < t>i^( r 2) ■ The amplitudes b[ are related to the spectroscopic factors by C 2 Si = \bi\ 2 . 
The wave function for > R is rather simple, it must be borne in mind that the binding itself is 
caused by the full three-body dynamics. 

The dipole operator is a one-body operator. It is given in terms of the Jacobi coordinates p\ and p 2 

by 

M{Elp) = Z cS e(p 1 Y lfl (p 1 ) + p 2 Y lfl {p 2 )) (4) 

The introduction of the Jacobi coordinates introduces an asymmetry between the two identical neutrons. 
This translates into an asymmetry of the effective charges for both neutrons, we neglect this effect, which 
is of the order of 1/A C , where A c is the mass of the core. For the effective charge we take Z e s = ~ z ^ l n . 
(The other possible choice would be Z e g = 2 ~m C ™m •) ^he final state is given by two continuum neutrons 

with momenta k\ and k 2 . We assume that they are uncorrelated. Low lying strength is a long wavelength 
phenomenon, and the interaction can be characterized by the n-core scattering lengths ao and a\ in the 
/ = and I = 1 channels. We know that 10 Li is unbound, so these scattering lengths are negative. The 
continuum radial wave functions are given by a product of single particle wave functions of the type 
given by eq. 3 above. 

The dipole matrix element consists of a product of two terms: one is characterized by the dipole 
selection rule Al = 1 and is of the type encountered in the single particle case above. We have 

dB(El) 2m 3 2 Cf lf 

where the shape functions [S] are given in terms of the dimensionless quantity x = ^ by 

^(a = i) = jyt^y^ 1 - a ^ + 3x2 ) + •■■) ( 6 ) 
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for s to p-transitions and by 



for p-to s-transitions. 

The other term is an overlap term 

dN 2 m , 

where E = 4^ (using the normalization conventions of [7]). One can describe the physics of this 
matrixelement in quite simple terms: Due to the sudden removal of one of the neutrons by the interaction 
of the photon with the core, the remaining neutron finds itself in the field of the 9 Li core, where it is no 
longer bound. A 'shake off' occurs , which is described by the overlap term, eq.9. It is non-zero, because 
the two wave functions belong to different Hamiltonians. A similar and familiar example is the shake-off 
of electrons after /3-decay, where the nuclear charge is suddenly changed by one unit. This approach 
can also be applied to non-Borromean nuclei. In this case the neutron can also be left in a bound state 
with a certain probability which can be obtained from the corresponding overlap matrixelement. 
The overlap integral is strongly peaked around zero energy. For I = one obtains 

< *o.l* 0t > I 2 = „ y fr*", 1 ?' (9) 



;i + M0 2 )(9 2 + k 2 )'- 



In the unitary limit (ao — > oo) this reduces to 



2ak 2 

< *o.|4V- > I 2 = (10) 

This is probably a good approximation for the overlap term in the n Li case where one neutron is 
almost bound to the 9 Li core. The analytical formula for all values of the angular momentum is given 
in Appendix A.4 of (A = 0). 

Low-lying El strength in n Li was observed experimentally in Ref. jHj. The dipole strength distri- 
bution as a function of the relative energy E = E 1 + E 2 = j-(kf + fcf) of the 9 Li-n-n-system is shown 
in Fig. 3 there. It will be interesting to compare this experimental result to the present approach. The 
dipole strength distribution E1 rf g n ~ ha '° depending on the relative energy E can be obtained by folding 
the dipole distribution with the overlap term: 

The factor of 2 arises because there are two neutrons in the same shell. We assume that the binding 
energy E b is known. The shapes of the functions Sq and S° (see eqs. 7 and 8) (and Sf) show charac- 
teristic differences, so the fit to the experimental results [H] will be quite sensitive to the spectroscopic 
factors C 2 S and C 2 Si. The width of the nonorthogonality term depends on the scattering lengths, so 
they can also be determined by this 4-parameter fit, but probably one will be less sensitive to these 
quantities. The experimental results [8] are quite well explained by the theory of Esbensen and Bertsch 
[T9*] . It will be interesting to compare the present approach also to the more detailed theory of [19J. 
In the present approach only a few low energy parameters like the scattering lengths in the relevant 
channels enter. 
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3 Trojan Horse Method 



A similarity between cross sections for two-body and closely related three-body reactions under certain 
kinematical conditions [2U] led to the introduction of the Trojan- Horse method [2H 1221 12H1 E] • in this 
indirect approach a two-body reaction 

A + x^C + c (12) 
that is relevant to nuclear astrophysics is replaced by a reaction 

A + a^ C + c + b (13) 

with three particles in the final state. One assumes that the Trojan horse a is composed predominantly 
of clusters x and b, i.e. a = (x + b). This reaction can be considered as a special case of a transfer 
reaction to the continuum. It is studied experimentally under quasi-free scattering conditions, i.e. when 
the momentum transfer to the spectator b is small. The method was primarily applied to the extraction 
of the low-energy cross section of reaction (|T2*|) that is relevant for astrophysics. However, the method 
can also be applied to the study of single-particle states in exotic nuclei around the particle threshold. 

3.1 Continuous Transition from positive to negative energies in the (A+x)- 
channel 

Motivated by this we look again at the relation between transfer to bound and unbound states. We 
study the reaction A + a — > B + b where a = (6 + x) and B denotes the final B = (A + x) system. It can 
be a bound state B with binding energy E hinA = —Ea x (> 0), the open channel A + x, with E& x > 0, or 
another channel C + c of the system B = (A + x). In particular, the reaction x + A — > C + c can have 
a positive Q value and the energy Ea x can be negative as well as positive. As an example we quote the 
recently studied Trojan horse reaction <i+ 6 Li [21] applied to the 6 Li(p, a) 3 He two-body reaction (the 
neutron being the spectator). In this case there are two charged particles in the initial state ( 6 Li+p). 
Another example with a neutral particle x would be 10 Be+<i — > p + n Be+7. The general question which 
we want to answer here is how the two regions Ea x > and Ea x < are related to each other. E.g., 
in Fig. 7 of j21| the coincidence yield is plotted as a function of the 6 Li-p relative energy. It is nonzero 
at zero relative energy. How does the theory j3] (and the experiment) continue to negative relative 
energies? With this method, subthreshold resonances can be investigated rather directly. 

The Trojan Horse amplitude f involves the combination / oc Si c ■ J ; + (see eq. (61) of Jlj). In this 
reference, only the case E = E^ x > is dealt with explicitly. The threshold behaviour of Jf in the case 
of neutrons is given by (A. 30), for charged particles by eq. (59) of [I], respectively. Si c is the S- matrix- 
element connecting the channel c = c + C and the elastic channel A + x with angular momentum 
For neutral particles the threshold behaviour is \Si c \ 2 ~ (kR) 2l+1 and the cross section 
tends to a finite constant. A similar analysis can be done for charged particles. Let us now study 
the case E& x < and the transition from E^ x > to E^ x < 0. In the latter case the wave number 
k becomes complex, k= iq and the factor involves a decaying function, instead of the oscillatory 
one. The quantity Si c is no longer an S-matrixelement, it is the normalization factor for the decaying 
wave-function. The continuous transition for both cases was studied in in a two-channel model. 

An interesting case is when there are subthreshold resonances, like in the systems lA N(p, 7) 15 and 
in 12 C(a,7) 16 0. In (21312] a two-channel model is studied. The standard Breit-Wigner result is obtained 

Sij = e^+^(<% V 3 . ) (14) 

3 v 3 E-E R + iT/2 J v ; 

where Er denotes the resonance energy, the partial widths are given by T\ and T2, the total width is 
given by T = Ti + T 2 . 
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For E < the width r\ turns out to be imaginary 9\. For En > we have a resonance, for Er < 
a subthreshold resonance, the formulae are valid for both cases. In l2Zj sub-Coulomb a transfer 
is used to study the low energy 12 C(a,7) 16 S-factor. Transitions to the 1~ and 2 + bound states are 
observed and the reduced a widths of these states is obtained. 



4 Conclusion 

While the foundations of direct reaction theory have been laid several decades ago, the new possibilites 
which have opened up with the rare isotope beams are an invitation to revisit this field. The general 
frame is set by nonrelativistic many-body quantum scattering theory, however, the increasing level of 
precision demands a good understanding of relativistic effects notably in intermediate-energy Coulomb 
excitation. 

The properties of halo nuclei depend very sensitively on the binding energy and despite the ever 
increasing precision of microscopic approaches using realistic NN forces it will not be possible, say, 
to predict the binding energies of nuclei to a level of about 100 keV. Thus halo nuclei ask for new 
approaches in terms of some effective low-energy constants. Such a treatment was provided in Ch. 2 
and an example with the one-neutron halo nucleus 11 Be was given. With the future radioactive beam 
facilities at RIKEN and GSI one will be able to study also neutron halo nuclei for intermediate masses 
in the years to come. This is expected to be relevant also for the astrophysical r-process. It is a great 
challenge to extend the present approach for one-nucleon halo nuclei to more complicated cases, like 
two-neutron halo nuclei. In the present work we obtained a rather simple formula for two-neutron halo 
nuclei with a transparent interpretation: after the first neutron is removed, the second one is left in 
a different field (Hamiltonian) which leads to a shake-off process, similar to the rearrangement of the 
atomic cloud after an a- or (3 decay. 

While the Coulomb dissociation method relies essentially only on QED, precise experiments can 
give, in combination with a thorough theoretical analysis, precise answers for the astrophysical S fac- 
tors. In the Trojan- Horse Method more phenomeno logical aspects enter, like optical model parameters 
and effective nuclear interactions. This makes the interpretation of the experimental results in terms 
of astrophysical S factors less precise. In the Trojan- Horse Method one can extract bare nucleus (un- 
screened) S-factors, see e.g. [21]. In the interpretation of screening effects one relies on the accuracy of 
the energy dependence. This is certainly better fulfilled than the accuracy of the absolute values. 
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